Bose-Einstein condensation in harmonic double wells 
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Abstract 

O 

T 7^ ' We discuss Bose-Einstein condensation in harmonic traps where the con- 

' finement has undergone a splitting along one direction. We mostly consider 



the 3D potentials consisting of two cylindrical wells separated a distance 2a 
along the z-axis. For ideal gases, the thermodynamics of the confined bosons 
has been investigated performing exact numerical summations to describe the 
major details of the transition and comparing the results with the semiclas- 
sical density-of-states approximation. We find that for large particle number 
and increasing well separation, the condensation temperature evolves from 
^ ■ the thermodynamic limit value T^\n) to (N/2). The effects of adding a 



repulsive interaction between atoms has been examined resorting to the Gross- 
Pitaevskii-Popov procedure and it is found that the shift of the condensation 
temperature exhibits different signs according to the separation between wells. 
In particular, for sufficiently large splitting, the trend opposes the well known 
results for harmonic traps, since the critical temperature appears to increase 
with growing repulsion strength. 

Pacs 03.75.Fi,05.30.Jp,32.80.Pj 
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I. INTRODUCTION 



The recent observation of Bose-Einstein condensation (BEC) of atomic alkali gases sub- 
jected to magnetic or magnetooptic traps (l|-|3| has triggered an important amount of related 
work. In fact, since the possibility of experimental realization of BEC was put forward, size- 
able theoretical effort addressed both the quantum and the thermodynamical aspects of 
this phenomenon when confined boson systems are involved. Starting from the earliest 
approaches to this field |j]-|7|], extensive work has been devoted to the various properties 
of harmonic traps [^|-|TU],|T^|-|2^] , which can account for most of the characteristic scales of 
experimentally achieved BEC according to very intuitive scaling arguments based on en- 
ergy balance [§] . The main features of BEC of ideal gases in either isotropic or anisotropic 
quadratic potential wells have been analyzed by various authors P Jl^JT2] -[T9| , paying special 
attention to the modifications of the transition patterns due to the finite number of trapped 
atoms. 

The nature of the transition itself has been the subject of some controversy. Kirsten and 
Toms have claimed that no spontaneous symmetry breaking leading to BEC can happen in 
confined potentials @,nj, as evidenced by the fact that the chemical potential of the bosons 
does not reach the value of the ground state energy at a finite temperature. However, 
this behavior of the chemical potential does not prevent the existence of a noticeable peak 
in the specific heat of the gas at a well defined temperature, for dimensions other than 
unity, coinciding with a sudden important rise in the ground state occupation number. This 
characteristic has been recognized by various authors [|lOH l2| as the signature of a phase 



transition. Furthermore, calculation methods based on exact summations 12-151 and Euler- 



MacLaurin approximation formulae |lqjr7|l yield results practically identical to those arising 
from more sophisticated treatments P, |10| ) p!8| . Bose condensed gases in anisotropic traps 
have also deserved attention, especially in view of the fact that experimental devices lead 
to the possibility that the definite stages of condensation are governed by the most binding 
ID forces |r7|,[nj. In particular, it has been clearly established that in highly anisotropic 
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potentials, the peak in the specific heat is linked to the freezing of those degrees of freedom 
lying higher in the energy spectrum, and macroscopic occupation of the ground state of the 
system occurs at a lower temperature . 



The effects of two particle interactions represented by a single parameter, the scattering 
length, have been examined P,|2"U|-f2"3]] ; the solutions of the Gross-Pitaevskii (GP) equation 
provide information upon the condensate wave function, and it has been shown that the 
transition temperature is sensitive to the sign of the interparticle forces in a way that opposes 
the expected results for the free homogeneous gas f22[. The variational solution to the GP 
equation p3| also provides an interesting means to approach the various features of BEC of 
weakly interacting gases. 

More recently, double Bose condensates have become an attractive field of research in 
view of the close resemblance between these systems and other bistable devices familiar in 



quantum optics [^4]]. The properties of double traps have been mostly investigated from the 
perspective of these analogies, keeping in mind that due to the large separations between the 
experimentally split traps and the significant height of the halving barriers, the condensate 
wave function essentially corresponds to two nonoverlapping wells. It is our present purpose 
to perform a detailed study of the quantum and thermodynamic aspects of these bistable 
traps in terms of the shape of the barrier and the well separation in the framework of 
a simple description. For this sake, in Sec. 2 we present the essential formulae for the 
quantum mechanics and the thermodynamics of an ideal gas in a harmonic double well, 
while the results of the calculation are presented in Sec. 3. Section 4 is devoted to the 
weakly interacting gas described by means of the GP + Popov formalism. The conclusions 
are summarized in Sec. 4. 

II. THE HARMONIC DOUBLE WELL 

The spectrum of the ID harmonic double well with potential 

V(z) = ^f(\ Z \-af (2.1) 
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can be represented [^5] by the parabolic cylinder functions D u [^2muj /h(\z\ — a)] with 
eigenvalues e v = fiu{y + 1/2)- The quantum numbers v are determined establishing that 
the eigenfunctions be either even or odd and it is found that the spectrum evolves from a 
purely harmonic one at separation a = to a doubly degenerate harmonic one of the same 
frequency , as a grows indefinitely. Furthermore, we shall consider 2D (i.e., uj v = 0) and 3D 
traps of the form 

V(r) = ^ [ulx* + uly* + ul(\z\-af] (2.2) 

A numerical computation of the energy levels of (|2.1| ) enables us to examine the thermo- 
dynamic properties of a confined boson gas by straightforward summation of the relation 
defining the internal energy, 

U = J £N i e i (2.3) 



where iV,- 



j3{Ei-fJ,) _ ^ 



is the boson occupation number at temperature T = 1/(3 for 



the level with energy E{ = Ti (n x u x +n y Uy + uu z + 3/2), once the particle number equation 

N = ]T n, (2.4) 

i 

has been solved to determine the chemical potential \x for the given value of N. The specific 
heat is afterwards computed differentiating the energy with respect to temperature. 

In this work we compare the results of the exact numerical calculation with those of 
the semiclassical method derived by Bagnato et al 0, which in the present case consists of 
replacing the summation in (|2.4j) by 

r _J e 

N = N + N 1 + Y,z J dep(e)e^T (2.5) 
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where p(e) is the classical extension of the density of states for single particles in the well. 
Notice that in Eq. ( |2.5| ) we explicitely separate the populations iVo of the ground state and 
Ni of the first excited state, since as the height of the barrier Vq = mu 2 a 2 /2 increases, these 
two energy levels approach the doubly degenerate ground state of the infinitely distant wells. 



The density of states is computed according to the prescription given in Ref. ||, which for 
an arbitrary confining D- dimensional potential V(r) can be cast into the most general form 

PD (e) = ^ 2 D ^ m D " J d D r [e - V^f^ 1 (2.6) 

where the integral is to be carried within the classically allowed region. In view of the 
relation 

PD {e) = f de' (e-eO Pl {e') (2.7) 
J o 

we only need to compute the ID density for the double well. The 3D case is studied 
convoluting p\ with 

P2(e) = ^- + (2.8) 

for uj = u x = uj y , which yields a formula for the number of particles in two dimensions that 
coincides with the high temperature limit of the exact summation (see Appendix). It should 
be noticed that the semiclassical prescription of Ref. |5| only gives the first term of Eq. ( |2.8| ) 
for the 2D density. 



III. THE SYMMETRIC DOUBLE WELL 

We have investigated BEC of an ideal boson gas confined in the trap described in the 
preceding section for dimensions D=l to 3, comparing the results of the exact summations 
with those of the approximate density-of-states method, in terms of the well splitting a 
and the number of particles N. In the 3D case, we also examine the anisotropic potentials 
mostly related to experimental situations; notice that in spite of the symmetry breaking 
along the z-direction, we shall consider the potential to be isotropic when 
and anisotropic otherwise. The analytical expressions for the semiclassical densities of states 
computed according to the prescription Q2.6p and the corresponding numbers of particles are 
summarized in the Appendix. 
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The specific calculations indicate that the ID well exhibits the characteristic sharp in- 
crease of the occupation number N around a temperature T c , however this is the usual 
seudotransition that is not accompanied by a peak in the specific heat, as shown in Fig. 
1 for various separations between harmonic wells (hereafter indicated by the dimensionless 



variable Zq = J2mco/fia). We can appreciate the moderate sensitivity of this quantity to 
the splitting z . We also compare the exact ground state ip( z ) of the double well with the 
wave function 

tp(z) = J\f [ip (z - a) + ip (z + a)] (3.1) 

where ipo(z ± a) represents the ground state of a single harmonic well of the same frequency 
centered at ±a and Af is a normalization factor. This is illustrated in Fig. 2, where we plot, 
for Zq =1 and 3, (left and right columns, respectively) the functions ip(z) and ip(z) in full 
and dashed lines, respectively. In the upper plots, each single ground state wave function 
ipo(z ± a) has been individually normalized, whereas in the lower ones the sum ( |3.1| ) has 
been normalized to unity. The dotted lines indicate each separate contribution to ip(z). It 
is clear that for every well separation here considered, the exact wave function resembles 
almost exactly the normalized combination (|3.1|) shown in the lower pictures. 



Let us now consider a 3D isotropic trap; no substantial qualitative differences have been 
observed in the 2D case. In Fig. 3 we exhibit the occupation numbers N and Ni as functions 
of the temperature in units of the critical temperature T^P' corresponding to infinite number 
of atoms, for well separations zq = (i.e., a single harmonic well) and 1, as well as for total 
particle numbers N=100 and 10000. We observe a rather sharp transition at a temperature 
slightly lower than unity for sufficiently large number of particles; furthermore, this transition 
temperature decreases with growing barrier height. On the other hand, the larger the value of 
Zq, the closer the resemblance between N and N\ even for temperatures substantially below 
Tj°); this fact reflects the quasidegeneracy of the ground state and the first excited state 
for large finite values of Zq, which evolves into complete degeneracy at infinite separation. 
This can be visualized in Fig. 4, where we have selected N=1000 and various barrier heights 
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Vo, which allow us to appreciate the significant decrease of the transition temperature that 
occurs as the wells split apart. 

The existence of a transition becomes evident as one analyzes the specific heat. Figs. 5 
and 6 correspond to the same parameters as Figs. 3 and 4, respectively; while lines indicate 
numerical results, symbols locate the calculations performed with the semiclassical expres- 
sions listed in the Appendix. We observe that the agreement between exact and approxi- 
mate results is almost perfect at the lowest temperatures, whereas for high temperatures, 
that agreement holds only for large particle numbers. For such numbers the semiclassical 
approach locates the transition to an acceptable accuracy, except for some sligth overes- 
timation of both the characteristic temperature and the height of the peak; however, this 
departure is very noticeable for small quantities of atoms. In Fig. 6 we may notice the 
consequences of increasing the well separation: the major effect takes place for a splitting 
Zo = 3, with only minor differences as one enlarges this number. 

We have also performed calculations for anisotropic traps; as an illustration, in Fig. 7 we 
show the specific heat for 1000 atoms and well separations Zq = and 1. The various curves 
correspond to the aspect ratios of the JILA [T|] and MIT traps 0, namely u z /u X: y = 
and {uj z /uj x = 3.2, uj z /uj y = 1.8), respectively, and to an isotropic well with a frequency 
uj = 0.56 uj z equal to the MIT geometric mean. We realize that the larger the anisotropy, 
the smaller is the transition temperature and the lower is the height of the peak. This also 
happens when one switches from the single to the double well, however to a less significant 
extent as compared to the anisotropy effect. 



IV. THE WEAKLY INTERACTING GAS IN THE DOUBLE WELL 

The influence of weak to moderate interactions between trapped particles can be in- 
vestigated in a GP mean field approach generalized for finite temperatures. Within this 



formalism, the condensate wavefunction ip(r) is obtained from the nonlinear equation |27 



ho^(r) = fi^r) (4.1) 



with 

h ° = {'J^ + V ^ + 9 [ ™ c(r) + 2 ^ (r)] } (4 ' 2) 

where we have made the usual decomposition of the density into condensate and noncon- 
densate contributions, i.e., n(r) = n c (r) + n(r), and n c (r) = \ip(r)\ 2 . The terms involving 
the interaction strength g arise from a two-body seudopotential g5(r). In the s-wave ap- 
proximation, which is adequate for very dilute gases, one has g = 4%ah 2 /m, being a the 
scattering length. 

Within the Popov approximation of the Hartree-Fock-Bogoliubov (HFB) theory, the 
excitations of a dilute system of bosons at low and intermediate temperatures are given by 
the coupled eigenvalue equations 

£ Ui(r) - g n c (r) Vi(r) = E { Ui{r) 

t Vi(r) - g n c (r) m(r) = -E { v^r) (4.3) 
which define the quasiparticle energies Ei and amplitudes Ui and Vi, with 

jC = h + gn c (r) (4.4) 



The noncondensate density is related to the quasiparticle properties according to 127,28 



n[r) 



) = E {K(r)| 2 + [k«| 2 + k(r)| 2 ] N{Ej] (4.5) 



where N(Ei) = (e@ El — is the Bose-Einstein distribution of the quasiparticles. 

One must solve the coupled GP (|4.1| ) and HFB (|4.3|) equations, using the self-consistent 
densities n c (r) and n(r) for fixed number of particles N. For this sake, we use the decoupling 
procedure introduced by Hutchinson et al |29| , and expand ho in a basis of eigenfunctions 
of the 2D harmonic oscillator on the (x, y) plane times a discretized basis set for the ID 
double well along the z-axis. To illustrate this procedure, in Fig. 8 we plot the condensate 
fraction as a function of the temperature, for 100 particles at various barrier heights and 



positive interaction strengths s = a Jmuj z /h. We observe that in all cases, the system 



reaches BEC at temperatures below the transition temperature for a single well with no 
interactions. However, depending on the value of z$ two different behaviors occur. When z 
is small, the condensation temperature decreases, the stronger the interaction, in agreement 
with the well known result for the harmonic trap |||2^]. Conversely, if zq is high enough, 
an increase of s yields a higher critical temperature. This feature can be interpreted if 
we consider the single particle spectrum in the mean field, whose ground and first excited 
states are depicted in Fig. 9 for a temperature T=0.26 and N=100. The overall trend is 
similar at any temperature and this choice is a convenient one in view of the displayed scale. 
We realize that these levels also exhibit two tendencies: for small well separations, as the 
interaction is enlarged the energy levels get closer, giving a more degenerate spectrum and 
a higher density of states, which in turn demand further cooling in order to achieve BEC. 
On the other hand, if Zq is sufficiently large (above roughly zq=2.5 in Fig. 9) the ground 
and first excited states of a double well are already quasi-degenerate; in this case, the effect 
of the interaction as indicated in this figure is to uniformly raise the whole spectrum. As a 
consequence, due to level crossing, particles occupying the first excited state of the free gas 
can be promoted to the interacting ground state at no thermodynamic cost, giving thus rise 
to a moderate increase of the condensation temperature. 



V. DISCUSSION AND SUMMARY 

In this work we have discussed the characteristics of Bose-Einstein condensation in har- 
monic wells where the confinement has undergone some splitting along one of the three 
directions. We mostly consider the 3D potentials consisting of two cylindrical wells sepa- 
rated a distance 2a along the z-axis. For ideal gases, the thermodynamics of the confined 
bosons has been investigated performing exact numerical summations for the particle num- 
ber, which enables us to determine the chemical potential, and for the total energy, after 
which we derive the specific heat by means of a numerical differentiation with respect to 
temperature. This procedure allows us to describe the major details of the transition and 
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to compare the results with the approximation corresponding to replacement of exact sum- 
mations by integrals weighted by a semiclassical density of states. 

The general results can be summarized as follows. We find that the semiclassical ap- 
proach is capable of locating the transition, i.e., the peak in the specific heat, with higher 
accuracy the larger the amount of trapped atoms. The overestimation in both the precise 
value of the transition temperature and the height of the peak is rather pronounced for 
small systems and becomes progressively less important as the system approaches the ther- 
modynamic limit. It becomes clear that the effect of halving the population is to lower the 
transition temperature; for large particle number, the trend of this temperature as the well 
separation grows from zero to infinity is clearly to sweep the path between critical tempera- 
tures T C (N) and T c (N/2). This general behavior is independent of the dimensionality of the 
system; however, just as in the single well problem fL2|]l4| , no strict phase transition takes 
place in the ID potential. In fact, in this case one finds that in spite of the important rise in 
the slope of the occupation number N observed at a rather well defined temperature, the 
specific heat increases monotonically towards the classical limit. 

The effects of adding a repulsive interaction between atoms has been examined resorting 
to the GP + Popov procedure which yields the condensate and noncondensate densities 
together with the single particle spectrum. It is found that the shift of the condensation 
temperature exhibits different signs according to the size of the separation between wells on 
the z axis; in particular, for sufficiently large splitting, the trend opposes the well known 
results for harmonic traps, since the critical temperature appears to increase with growing 
repulsion strength. This unexpected feature can be explained examining the spectrum of 
the atoms in the mean field as a function of the well separation. One can realize that in the 
interacting system, not only the position of the energy levels, but the size of the gap between 
the ground and the first excited state, change substantially when the wells split apart. The 
level crossing that occurs for large separations causes atoms in the excited state to move 
into the new ground state, with a consequent increase in the condensation temperature. 

An experimental realization of the double well is the MIT trap characterized by a 
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large separation between wells and by a tunable barrier height. One might then wonder to 
what an extent the results here presented might be sensitive to independent choices of the 
location ±a of the minima and the barrier magnitude Vo, and intend to develop more detailed 
models for the specific MIT double trap. Among the simplest choices, one could quote i) a 
combination of two symmetric harmonic wells smoothly joined at z = ±a m by a quadratic 
barrier and ii) the double well of the previous sections deformed by a central constant cutoff 
with value Vo for \z\ < a m . We have performed calculations similar to those presented 
in section 3 for the above potentials, employing standard magnitudes of the original MIT 
trap |§ . We have found that the thermodynamics of the ideal gas remains identical to the 
previously analized double well case, in other words, the transition temperature and the 
overall shape of the occupation numbers and specific heat (Figs. 3 to 6) are not sensitive to 
the decoupling between the location of the minima and the barrier size, within the range of 
values of interest. 

To summarize, we remark that the statistical mechanics of a bistable system can be 
examined on identical grounds as the case of a single equilibrium state, with an additional 
parameter, namely the separation between wells, as the agent of the evolution from a given 
number of trapped atoms up to the point where one exactly halves the population in the 
infinite separation limit. 
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APPENDIX 

We collect here various useful formulae concerning the thermodynamics of ideal gases in 
harmonic traps. We display here the semiclassical results obtained from the density-of-states 
method of Ref. |§ (Eq. (|2.5| ) for the isotropic D-dimensional oscillator, 
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Particle number: 



D 



Internal energy: 



fkT\ D+1 



Specific heat: 



C 



'kT 



D 



D . 

k \ Tiuj 



(D + l)g D+1 (z)-D^l(z) 

9D-1 



whereas the exact summations give, in the high temperature limit hoo << kT 



kT 



d D 



z e 



{n-D/2)hu/kT 



(5.1) 



(5.2) 



(5.3) 



(5.4) 



In this limit one can take advantage of the property of the Bose special functions gi(z) 

12 + 5/2 

<Ji\e kT I » //,(:) + 



Sfj, gi-i(z) 



kT z 



(5.5) 



for D > 2 and write 



to(^) + — T— 



D kT g D -i(z) 



(5.6) 



2 hu z 

from where the transition temperature is obtained specifying N = together with z — 1. 
For the isotropic potential with the double well on the z direction, we have the following 

cases 



Dimension D — 1: 

The density of states takes the form 



Pid{s) = < 



J2_ 



1,21 



(5.7) 



and the number of particles is 



12 



+ de 



Z + Zl + In _ Zm 
1 — 2 1 — z\ (3fao (1 — 21] 

2 1 r 00 /K 
H — - — / dearcsin\ — N(s) 

7T «,Co> J £m V £ 



(5.8) 



with z = exp (/?//), = exp {f3(fx — ei)} ,z m = exp — e m )} and e m = maxjei, K} The 
total energy can be expressed, similarly to (|5.8| ), as a summation including an integral. 



Dimension D=2: 
We obtain 

P2d{s) = 

and 

iV 



(hu) 



£ <V n 



e + ef arcsinUVo/e) + f yJV yJe - V Q , e>V Q 



Z Z\ 1 r , . , 

+ H "T7TZ To [2 Cfe(*l) - </2(Zn 



l-z l-z 1 ((3hu) 2 
2 1 



+- 



IT (3(kuJ 

where g n (z) is the usual Bose function [26] 



, de arcsin \ hi 1 — cxp( — ?(< — // ))| 



(5.9) 



(5.10) 



Dimension D=3: 

The corresponding expressions are 



and 



P3D(£) 



IT- 



£<V 



e 2 + ^e 2 arcsin ( J~^f- 



£>V 



(5.11) 
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N = + j^— + -^-^ {2 g 3 ( Zl ) - g 3 {z m ) + (3hu [2 g 2 ( Zl ) - g 2 (z m )]} 

i — z 1 — Z\ \jjaui)° 



+ 



2 1 f°° ( fv~\ 

KP^L d£arCSm [ii ^(Ze) + Phu 9l (z £ )] (5.12) 
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FIGURE CAPTIONS 

Figure 1.- The specific heat of a ID boson gas trapped in a double well (in units of the 
Boltzmann constant k B ) as a function of the reduced temperature T/T^ for different 
separation between potential minima. 

Figure 2.- The exact ground state probability [V^)) 2 (full lines) of particles in the double 
trap as a function of the dimensionless coordinate z/a, for well separations zo=l (left 
column) and z =3 (right column). Dashed lines correspond to the sum of two normal- 
ized single well wave functions (upper plots) and the normalized sum (lower plots), 
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while dotted lines indicate the individual ground states. 

Figure 3.- Occupation numbers of the ground and first excited state for total particle 
number N=100 (thin lines) and 10000 (thick lines) as functions of T/T^ for the 
single harmonic well and for a separation equal to unity. 

Figure 4.- Occupation numbers of the ground and first excited state for 1000 atoms and 
various barrier heights. 

Figure 5.- Specific heat for the same situation depicted in Fig. 3. Circles and dots indicate 
the results of the semiclassical approximation. 

Figure 6.- Same as Fig. 5 for the conditions of Fig. 4. 

Figure 7.- Specific heat of anisotropic traps for the aspect ratios of the JILA and MIT 
traps, compared to an isotropic situation. 

Figure 8.- Ground state occupation number of a 3D double well including weak two-body 
interactions, computed in the GP+Popov approximation. 

Figure 9.- Energies of the ground and first excited states (full and dotted lines, respec- 
tively) of interacting particles in an isotropic double well, as functions of the separation 
Zo, for various interaction strengths from s=0 up to s=10~ 2 (bottom to top) in steps 
As=10~ 3 . The dashed line indicates the barrier height V . Energies are given in units 

of huj. 
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